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Abstract. We consider three-body systems in two dimensions with zero-range 
interactions for general masses and interaction strengths. The momentum- 
space Schrodinger equation is solved numerically and in the Born-Oppenheimer 
(BO) approximation. The BO expression is derived using separable potentials 
and yields a concise adiabatic potential between the two heavy particles. The 
BO potential is Coulomb-like and exponentially decreasing at small and large 
distances, respectively. While we find similar qualitative features to previous 
studies, we find important quantitative differences. Our results demonstrate that 
mass-imbalanced systems that are accessible in the field of ultracold atomic gases 
can have a rich three-body bound state spectrum in two dimensional geometries. 
Small light-heavy mass ratios increase the number of bound states. For 87 Rb- 
87 Rb- 6 Li and 133 Cs- 133 Cs- 6 Li we find respectively 3 and 4 bound states. 
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1. Introduction 

Manipulating quantum systems on large, small and intermediate scales is no longer 
a distant dream. Ultracold quantum gases experiments have realized Bose-Einstein 
condensates pQ and degenerate Fermi gases [5] and study truly remarkable many-body 
states, but it has also proven its worth in exploring the physics of few-body systems 
[3]. Moreover, these studies can be conducted in different geometries as two- and 
one-dimensional quantum systems are produced regularly in laboratories pQ. 

The so-called Efimov three-body bound states were predicted over four decades 
ago [1]. They arise when a three-body system, composed of equal mass particles, 
has all its two-body subsystems at the threshold of binding. In this case, infinitely 
many bound states with energies geometrically separated are expected. These are 
the Efimov states and this effect is called the Efimov effect. The effect is related 
to a long-range effective force and it can occur even when the individual two-body 
forces have zero range, as it was shown in an analytically solvable model by Fonseca, 
Redish and Shanley J5, . This is an example of how long-range forces can arise in the 
three-body problem in a way unpredictable by two-body intuition. In the early work 
of Ref. [5] the Efimov problem was handled through the adiabatic approximation, 
more precisely the Born-Oppenheimer (BO) approximation that is commonly used in 
molecular physics. The BO approximation considers a system composed of two heavy 
and one light particle. The terms heavy and light have relative meaning: two particles 
are heavier than the third one. In this approximation the heavy particles move very 
slowly while the light particle orbits around them. In fact, it is enough to assume 
that the heavy particles kinetic energy is (much) smaller than the kinetic energy of 
the light particle. 

These results relating the Efimov effect and the adiabatic approximation were 
obtained for three-dimensional (3D) systems. The dynamics and properties of 
quantum systems drastically change when the system is restricted to lower dimensions. 
For example, the scattering-length is not well defined for two-dimensional (2D) systems 
[6] and the kinetic energy operator gives a negative centrifugal barrier for 2D systems 
with zero total angular momentum while the centrifugal barrier is always non-negative 
for 3D systems. Furthermore, it was shown that any infinitesimal amount of attraction 
produces a bound state in 2D [TJ [HI El US], while a finite amount of attraction is 
necessary for binding a 3D system. 

Another important difference between 2D and 3D systems is the occurrence of 
the Efimov effect. The effect does not occur in 2D neither for equal mass three- 
body system [TT| or unequal mass systems [IS]. While the BO approximation 
was implemented in Ref. [IS] to look for the Efimov effect in 2D imbalanced mass 
three-body systems, the mass-dependence of such systems was not addressed. The 
importance of the mass dependence in the 2D imbalanced mass three-body systems 
was stressed in Refs. [El [18], where an increasing number of bound states was found 
for a decreasing mass of one particles. This situation, where one particle is much 
lighter than the other ones, is suitably handled in the adiabatic approximation. 

In this work we study the BO approximation of 2D three-body systems from the 
mass-dependence perspective in a systematic fashion. We consider a 2D three-body 
system with zero-range interactions for general masses and interaction strengths. The 
momentum-space Schrodinger equation is solved numerically both in full generality 
and in the BO approximation. The BO expression is derived using separable potentials 
and yields a concise adiabatic potential between the two heavy particles in the heavy- 
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heavy-light system when the light particle coordinate is integrated out. The adiabatic 
potential is mass-dependent and it reveals an increasing number of bound states for 
the decreasing mass of the light particle. We find a transcendental equation that 
describes the adiabatic potential and an approximate analytic expression as well as 
the asymptotic form. Furthermore, we show that the approximate analytic form of 
the potential is very close to the full adiabatic potential. 

Furthermore, we estimate the number of bound states for a heavy-heavy-light 
system as a function of the light-heavy mass ratio. Infinitely many bound states are 
expected as this ratio approaches zero. However, for each given mass configuration, 
we still have a finite number of bound states. Our results demonstrate that mass- 
imbalanced systems that are accessible in recent ultracold atomic gases experiments 
can have a rich three-body bound state spectrum in 2D geometries, where small mass 
ratios increase the number of bound states. For 87 Rb- 87 Rb- 6 Li and 133 Cs- 133 Cs- 6 Li 
we find respectively 3 and 4 bound states. 

The paper is structured as follows. The introduction in SectionfTJis followed by the 
2D three-body system formalism for zero-range interactions and general masses in the 
BO approximation in Section [2] The BO approximation and the adiabatic potential 
are shown in Section |3l The results are presented in Section |4] and a discussion is 
given in Section [5] 

2. Formalism 

The general setup that we consider has particles that are always confined to a two- 
dimensional plane of motion, i.e. the kinematics is always 2D. However, we do not 
rule out systems with long-range interactions across many two-dimensional planes (as 
done in recent experiments with polar molecules [13]). Both of these situations are 
shown in figure Q] for the most general case with three different particles A, B and 
C . The only assumptions concerning the interactions is that they are dependent only 
on the relative distance of the particles under consideration. We can then assume 
that the three-particle dynamics effectively happens in a single plane, since we may 
consider the layer index as an additional quantum number that specifies a few-body 
state. 

Here we are considering only zero-range interactions since we are interested in 
the model-independent low-energy universal limit. This provides a simplification 
of the Faddeev equation for the three-body bound state due to the separability of 
the zero-range interaction. In the case where we have long-range interaction across 
different layers and we have to assume that at low energy, the long-range forces 
may be approximated by a short-range pseudopotential that reproduces the two- 
body properties of the full potential. As discussed in Ref. [14] . this can be done 
for potential that decay faster than r~ 2 for r — > oo. This includes the case of dipolar 
r~ 3 interactions that we will return at the end of this paper. 

The 2D Hamiltonian for the three-particle ABC system with pairwise two-body 
potentials is 

H = H Q + V AB + V AC + V BC (1) 

where the three particles are, in principle, assumed to be bosons. However, for three 
different mass particles there is no symmetry requirements. Our formalism takes the 
symmetry fully into account when two particles are identical bosons. The kinetic 
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Figure 1. A pictorial representation of the setup under consideration, (a) A 
three-body system in a plane (viewed from the top). We show the plane as 
a pancake with ellipsoidal shape rather than as a square depiction in order to 
convey the fact that real system may have deformation of the trap that keep the 
particles fixed in two dimensions. The three distinct particles are labeled A, B 
and C. The two-body binding energies, Ea_b, Eac and Ebc> and the coordinate 
axis are shown, (b) A multi-layer setup with several A,B and C particles (side 
view). Long-range interactions in and between the layers can lead to bound state 
formation both within the same layer and for particles in different layers. In the 
case of three-body states, this is indicated by black arrows. 



energy operator is 



PA 



q| 



p! 



q 2 c 



P 2 c 



2/M,BC %fJ-B,AC 2fJ-AB ^^C,AB 2fJ-AB ' 



(2) 



where we use Jacobi relative momenta given in terms of rest frame momenta, with 

i = a, ft, 7, as 

q 7 = k 7 and p 7 = fj, a p \ — — 

\m a mp t 



(3) 



where (a, /3,7) is the cyclic permutations of the particles A, B, and C with masses 



m,A, mB and mc- The reduced masses are ix a p — 
Wc use a separable potential with operator form 



and /i 7 , a ^ 



V a p = X a /3\Xa0}(Xap\, 



(4) 



where the form factor is (p 7 | Xap) = gip-y), depends only on the relative momentum 
of the two particles a and /3. The limit to zero-range interaction is simply g(p) = 1. 

In order to relate the separable potential ansatz above to physical properties of 
the two-body system we use the the two-body T-matrix. For negative energies and 
zero-range potentials it is defined by 



T a p{E) = \x a p)T a p(E)( Xa p\, 



(5) 
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where, the matrix element of the 2D transition matrices are given by (see e.g. Ref. |15j 
for the case of identical particles) 

-l 




(6) 



where a,j3 = A, B or C and a ^ (3. E a p is the energy of the a(3 two-body bound 
states. We adopt units such that h= 1. Also we will measure all three-body energies 
in units of the two-body energy. In general systems where A, B, and C are different, 
one only needs to make a specific choice between the E a p. We will explicitly state 
our units in these cases. Also it is important to know that E a p = it corresponds 
to a non-interacting a/3 system since T a p — > in this case. This is different from 
3D where a two-body bound state at zero energy happens at unitarity, i.e. where 
the interactions are as strong as allowed by the unitary limit of quantum mechanics 
applied to the s-wave scattering amplitude. 

In order to solve for three-body bound states, the wave function \^abc) is 
decomposed in terms of the Faddeev components as 

\*ABc) = + |*B> + l*C> (7) 

where |* 7 ) = G Q (E)V a p\^> abc) with the resolvent G (E) = (E - Hq)^ 1 which is 
nonsingular for bound states. The bound state equation can be solved by the method 
described in Ref. [TB] for AAB systems with some straightforward generalizations to 
the ABC structure. We therefore refer to Ref. [TB] for the details of the formalism 
and do not repeat them here. Instead we now proceed to discuss the adiabatic 
approximation . 



3. The Adiabatic Approximation 

We now consider two heavy particles with mass vtla and tub- These particles are 
fixed and their centers are separated by a distance R. The additional light particle 
has mass mc and coordinate r relative to the center-of-mass of the heavy subsystem. 
The configuration and coordinates are depicted in figure HJ We assume that the 
particles interact with each other through short-range potentials. The notation for 
the potential is that vc means the interaction between particles A and B and va and 
Vb are analogously defined. The three-body Hamiltonian is given by 

2 2 

H = + - q °' +va+vb + vc , (8) 

l[lAB t\lC,AB 

where the Jacobi relative momenta and reduced mass definitions are given in ([3]). The 
Schrddinger eigenvalue equation i?^(r, R) = E^fr, R) where is given by 



ff = -^V|--^V? + ^(r-^R)+, B (r+^R]+ Uc (R)- 

2fj, AB *nc,AB \ m B J \ m A 



(9) 



The adiabatic approximation now instructs a split of the three-body eigenvalue 
equation into the solution of two two-body problems: the light particle motion is 
considered with respect to the heavy-heavy system and the heavy-heavy system motion 
is separated out. These eigenvalue equations will be useful whenever the motion of 
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Figure 2. Three-body relative coordinates used in the adiabatic approximation. 
Heavy particle of mass vixa and ms interact with a light particle of mass vac- 
The heavy particles are fixed at a distance R from each other and the origin of 
the coordinate system is chosen as the center-of-mass of the two heavy particles. 
The coordinate of the light particle is its distance from this choice of origin. 



the light particle is rapid compared to the motion of the heavy ones, so that the light 
particle dynamics can be solved while the heavy particles are instantaneously at rest. 
The wave function is written as 



tf(r,R)=^(r,RMR) 



(10) 



where ip(r, R) is the wave function describing the motion of the light particle for fixed 
R and 4>(R) is the heavy system wave function. The approximation is valid when the 
kinetic energy, — 2 V^V( r i R); is small compared to the other terms in ((9|). Using 
the wave function (|10p . we write the light particle eigenvalue equation as 

h 2 _ 



' V A 



AB 



HAB 

m B 



R 



VB 



fJ-AB 

m A 



R 



V>(r,R) = e(i?)V>(r,R). (11) 



where e(R) is a separation constant, which according to this equation cannot depend on 
r. Furthermore, e(R) plays the role of an effective potential in the equation describing 
the heavy-heavy particle system, i.e. 



2fJ,AB 



vc{R) + e(R) (f>(R) = E<j)(R) 



Assuming that the potentials in (|lip are separable and have the same strength, i.e. 
v a = A \xa) (Xa\-> the light particle wave function in (jlip in momentum space reads 



where 
A+ 



<l>(p) = X 



d 2 r'~g^ 
d 2 r'g^ 



e(R) 



^ p 2 



-4, 



I-IAB 

m A 

(J-AB 

m B 



RJ * (r') 
R ] * (r') 



<W (P') 



<W (p') 



^p'R 



2tt 



-p' R 



2tt 



-^(p') 



(13) 

(14) 
(15) 
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Rewriting (| 13[) in terms of A± , we find the system of equations 



A± = A / d 2 P 



is(p)r 



e(R) 



(e Tip - R A T +A ± ) 



(16) 



The non-trivial solution of (JT6J) , i.e A± 7^ 0, gives a transcendental equation for the 
energy, which reads 



1 
A 



d 2 p 



I.9(P)| 2 



e(R) 



h 2 



[1 + cos(p • R)] 



(17) 



We now use the binding energy of the two-particle heavy- light subsystem, Ei , and the 
T-matrix to eliminate A [5] . This yields 



d 2 P \g(p)f 



1 + cos(p • R) 



^p 2 \E2\ + 



2fJ-C,AB 







(18) 



Since we are interested in low-energy, model-independent bound state, we describe 
the heavy-light particle system interacting through zero-range interactions. In 
momentum space this means that <?(p) = 1, so (|18|) is finite and its two terms can be 
integrated to give the transcendental equation 



'2hc,ab\<R)\ 



n 2 



R 



(19) 



where K$ is the zero order modified Bcsscl function of the second kind. The (119]) 
is a powerful tool in understanding mass-imbalanced three-body systems in two 
dimensions. We are able to solve it to get an analytic form for the small and large 
distance behavior of the adiabatic potential. We can also numerically solve it to test 
the validity of the analytical expressions. 

When the separation, R, between the two-heavy-particles is large, i.e. |R| — > 00, 
the light particle feels only the interaction from one of the heavy particles. In this limit 
the three-body problem becomes a two-body problem and we expect that \E\ = \E2\- 
Thus, defining Jf^i = 1 + V(R), this condition is fulfilled when V -> for |R| -> 00. 
Inserting these definitions in (|19l) and expanding both sides up to first order in V(R), 
we obtain 

2K (s(R)) 



V(R) 



l + s(R) Ki (s(R)Y 



(20) 



where s(R) — y 2 ^ C ?^I- E3 1 \ n f ac t ; ([20]) is a very good approximation of the BO 
potential e(R). We are also able to solve (fT9|) to get the small distance behavior of 
the adiabatic potential. Using the asymptotic form of Kq for small arguments |17) . 
the adiabatic potential reads 



\E 2 \ 

for R -> and 



2e~~t 



e~ 7 
1 - — s(R) 



(l- 7 ) --In 



-s(R) 



(21) 



(22) 
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Once the adiabatic potential is known, we can solve (|12[) and determine the bound 
state spectrum. In order to write a Sturm-Liouville eigenvalue equation for the heavy- 
heavy system in a L z = state, we set <f> = -j= and (fT2")l becomes 



with the adiabatic potential from (fT9")l given in the limits of R — > and R — >• oo by 
(|2Tj) and respectively. 

Due to the behavior of the adiabatic potential at small and large distance, we 
expect to find an increasing number of bound states when particle C is much lighter 
than the other ones (mc —> 0). In this limit the adiabatic potential becomes more 
attractive and these states will accumulate both at R — > as the strength of the 
Coulomb potential increases, see (f2"Tj) . and at R — > oo, where more states are allowed 
because the exponential moves to larger distances, see (|2"2"|) . However, for finite mc, 
the number of bound states is still finite. 

One might argue that the limit \E2 \ — > must produce the same effect as mc — > 
in the adiabatic potential of (f2~Tj) and (|22l) . The limit where all subsystems interacting 
through zero-range interactions are unbound does not support three-body bound states 
[TTJ [T51 [T5] . However, recent studies show that these so-called Borromean states are 
possible in 2D but only for potentials with an outer barrier and an inner attractive 
pocket pH[20]. 

4. Results 

Adiabatic potentials 

We found analytical forms for the small, (|2"Tj) . and large, (1221) . distance behavior of 
the adiabatic potential. Once we are able to numerically compute the exact adiabatic 
potential from (fl9|) , we must show that our asymptotic behaviors in (|2"Tj) and ([22| are 
consistent with the exact adiabatic potential. 

The exact potential, e e xact(s), obtained numerically from (|19[) is shown in figure [3] 
and compared to e a symp(s) of Ecis. [2Tl and (|22|) . The asymptotic behavior does indeed 
work very well for most of coordinate space. The largest deviations are found in the 
region 0.3 < s < 3, where the difference between e asym p(s) and e exact {s) still never 
exceeds 9%. Notice that the approximation accuracy increases a lot when higher 
order terms are included in the expansions. One could go to more precise adiabatic 
potential representations taken higher order expansions of (|19p . However, the results of 
the approximate and full adiabatic potentials are almost indistinguishable in practice. 

4-2. Bound states 

It is known that the number of bound states increases as one particle becomes much 
lighter than the other ones [TBI HE]. Here, we have calculated the adiabatic potential 
asymptotic forms (|2T]) and ([22]) and found a differential equation ((23]) applicable to 
the case where two particles are much heavier than the third. In the following, we 
specialize to identical heavy particles and adopt units such that E^c = Egc — ^2 
and m^ = mg = M = I. The mass ratio between light and heavy particles is then 




(23) 
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S -> 00 

s -> (1 sl order) 

■ - -s -> (2 nd order) 
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Figure 3. Ratio £ a symptotic(s)/ e exact(s) as function of the dimensionless 
coordinate s, showing the validity of the asymptotic behavior in (12 1 1 ) and {22) . 
The solid (black) and dashed (red) curves are the first order expansion at large 
and small distances, respectively. The dotted (green) curve is the second order 
expansion (12 U . 



denoted m = as in [16]. In this case, the reduced mass hc,ab is written as 



2m 



(J-C,ab 



and [ic.AB m for m — > 



(24) 



With these dehnitions and denoting the three-body energy by E = £3, a useful 
adiabatic potential expression is given by 



e(R) -> - 



2e"T 

4m td 
m+2 n 



1 - 



Am 
m + 2 



R 



1 



(1 - 7) - - In „ 

V " 2 I 2 V m + 2 



4m 



R 



(25) 



for ,/-^2- R < 1.15 and 

m+2 — 



e(R) -> -1 



27T 



for 



Am 
m+2 , 



R 



4m 



-R > 1.15. 



(26) 



This approximation is very accurate when 2R « 1.15^/(1 + 2/m) where the largest 
deviation of 9% is reached. Proper bound three-body states are present when 
E 3 - E 2 < 0, or equivalently \E 3 \ > \E 2 \. 



We numerically solve the differential equation (|23|) with the adiabatic potential 
(|25p and (J26J to estimate the number of bound states (Nb) for a system with mass 
ratio m when the heavy particles do not interact with each other. Letting Eab 
be the binding energy of the heavy-heavy system, we can choose vq = in (j23[) 
and get E^b = 0. If t>c is attractive and able to support bound states, the three- 
body system would effectively be reduced to the lightest particle moving around a 
heavy-heavy dimer. The corresponding additional much deeper-lying bound states are, 
however, not interesting in the present context. Such states are not universal three- 
body structures that we address in the current study, but should rather be regarded as 
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bound two-body states dressed by an orbiting third particle. The numerical solution is 
obtained by writing the eigenvalue equation (|2l?)) in matrix form. After discretisation 
of the operators and the wave function of the radial equation, we get a tridiagonal 
matrix which can be diagonalized to obtain the energies and the number of bound 
states. 



K>(> 




10' 1 10 J 10'' 10' 1 10'' 10° 

m 



Figure 4. Number of bound states, Nb, for a system with mass ratio m and 
Ej\g = 0. The (black) squares represent the mass ratio m where Ng states are 
bound, calculated from the adiabatic approximation in i'2'it . The (red) circles 
represent the full numerical solution. The inset shows the result for the largest m 
values on a different scale. 



The number of bound states as function of mass ratio, m, is shown in figure HI 
where a comparison between the adiabatic result of (|23|) and the full numerical solution 
is made. Comparing both results in figure [4j we see that the adiabatic approximation 
picks up the small mass behavior very well, even for mass ratios up towards 1. There is 
a small error in the threshold for the number of available bound states for < Ng < 14, 
but it decreases as m — ¥ 0. The adiabatic approximation has an accuracy better than 
10% for m = 0.01, as can be seen in the inset of figure |3J Due to the numerical 
difficulties, it is very hard to count the number of bound states for Nb > 14 by 
solving the full problem numerically. Fortunately, it is very easy to do it through the 
differential equation (|23|) with the adiabatic potentials ([25]) and (|26| . We also see 
clearly in figure |4] that Nb — > oo for m — > 0, as pointed out in [16]. 

A fit to the results presented in figured shows that the dependence of the number 
of bound states, Nb, with the mass ratio, m, is rather well described by 

N B » "S. (27) 



This behavior can be explained by the old quantum mechanics. What we need to 
consider is the number of nodes in the wave function at zero energy, since this measures 
the number of allowed bound states. Within the old semi-classical quantum theory, 
the usual way to estimate the number of bound states in an one-dimensional quantum 
problem is 



p dq = NttH. (28) 
Taking into account the effective potential in ()23[) and proper units, the number of 
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bound states can be estimated from the formula 

1 



N 



dx 



m 
2x* 



V(x) 



0.733 



(29) 



where V(x) is the adiabatic potential (|23|) and (|2"B1) with x = y ^^R- One may 

object that the integral in (j2"9")l diverges in both limits and can not be performed. 
Introducing a lower and an upper cut-off in the integral, which are the same used in 
the numerical calculation (1CP 2 and 10 5 respectively), we get that N — ■ This 
result approaches the estimate in (|2"T|) once the diverging term in (j2"9")l becomes less 
important as m becomes smaller. The integral in (f29| is m-independent for m < 0.001 
with a 10~ 2 cut-off, implying that the term ra/x 2 is negligibly small itself. The 
apparent divergences are due to the semi-classical estimate, and accurately removed 
by a cut-off at both small and large x. The true quantum mechanical number of states 
can then be fully recovered. 




10 10 10 10 10 10 10" 10 10 



Figure 5. Number of bound states, Ng, for a system with mass ratio m for 
Ej^g = (dashed (red) line) and for general (but fixed) Ej^g (full (black) line). 
The (black) squares represent some of the points from figure [4] 

The estimates of the number of bound states in ([27)) and (|29j) agree very nicely. 
These estimates are less than the upper two-dimensional limit for a two-dimensional 
system with total angular momentum equal to zero, which is given in Ref. 21j. For 
the adiabatic potential in ([23)1 and (|2"o| . this upper limit is given by N = — . The 
difference between the estimates is shown in figure [SJ It is known that any three-body 
system in two dimensions will achieve its maximum number of bound states when all 
two-body subsystem are bound and have the same two-body binding energy [18] . We 
thus expect that the estimate given by the dashed (red) curve in figure [5] will hold for 
the adiabatic potential in (|25|) and (|26|) when Eab — E2. Also, the number of bound 
states for a system with < I-E12I < 1 is in the window between the solid (black) 
curve and the dashed (red) curve shown in figure The large difference between 
these curves is necessarily due to states that have non-zero angular momentum that 
we are not concerned with in the present study but that are included in the strict 
mathematical bound cited above. 

As we expected, the results confirm that the bound states accumulate at both 
R — > and R — ¥ 00 as m — > 0. From our numerical calculations we find that the 
energy of the lower states (which reside at small distance) seem to increase without 
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bound when m — > 0. Furthermore, in the same limit the wave junction vanishes 
slower at large distances, allowing more bound states also at large distance. This can 
be interpreted as an Efimov-like effect for the two dimensional case. However, an 
important distinction between the 2D and 3D case must be made. While the Efimov 
effect says that three identical particles can have infinitely many bound states when 
E% — > 0, in 2D this limit leads to an unbound three-body system. We expected to 
have infinitely many bound states in 2D, when m = and two interactive interactions. 
But we must stress that for finite m we still have a finite number of bound states. 

4-3. Experimentally relevant systems 

Experiments with ultracold atoms are able to produce quasi-2D samples of 23 Na [22] , 
40 K [231 ED, 87 Rb [25], 133 Cs [HI EH [28], and 6 Li [59]. It was recently reported 
that mixtures of 133 Cs and 6 Li were successfully trapped and that very favourable 
Feshbach resonances that can be used for tuning the interaction strength between the 
atoms have now been found [3D1 OH]- In order to discuss some properties of mass- 
imbalanced three-body systems, we consider two different systems that can be probed 
in laboratories in the near future. These system are 133 Cs- 133 Cs- 6 Li, represented as 
circles in figures [7J and [5] and 87 Rb- 87 Rb- 6 Li, represented as diamonds in the same 
figures. 



1 2 3 




Figure 6. Pictorials of the relation between masses and energies for regions 1, 2 
and 3 of shown in figures [7] and [8] 

Figures [7J and [5] are mass-diagrams that show how many bound states a system 
composed of different particles supports for a given set of two-body binding energies 
in each subsystem. The energy symmetric case (all two-body binding energies equal) 
was discussed in Ref. [18]. In general there are several ways to obtain a 133 Cs- 133 Cs- 
6 Li system if we assume that the interactions in the different subsystems give rise to 
different binding energies. The possiblities are shown in figure [T5] If we consider a 
system of just one hyperfine state of 6 Li and of 133 Cs clearly two of the interactions 
will be identical (for instance Eac — Ebc on the left side of figure H]). However, as 
has been discussed in the recent studies of Feshbach resonances in this mixed system 
[30] [31] different hyperfine states will in general give different interactions. This implies 
that with different hyperfine states we may be able to access the full range of parameter 
space. The adiabatic problem was handled for a non-interacting heavy-heavy system, 
i.e. Eab — 0. This is close to the situation in 133 Cs- 6 Li experiments, where three- 
body bound states are expected to be found when the subsystem 133 Cs- 133 Cs is almost 
non- interacting [30l[3T]. The relevant case is then Eab = which is shown in figure [Jj 
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Figure 7. Mass-diagram defining the mass ratios thresholds for the existence 
of N bound states for Ej\g = 0,Eaq = EgQ. The numbers 1, 2 and 3 label 
three different regions according to figure \E\ and the Roman numerals denote the 
number of bound states present bellow each solid curve. The circles mark the 
133 Cs- 133 Cs- 6 Li systems and the diamonds mark the 87 Rb- 87 Rb- 6 Li systems. 



The mass-diagram in figures [7J and [5] was obtained by solving the full numerical 
integral equations as in Ref. [TB] . In figure the points at which the lines cross 

= 1 corresponds to the red circles in figure 2] Notice that on region 1 , where 
mc < mA,mB, the particle C is always the lightest and Eab is always the heavy- 
heavy particle energy. This region gives the deepest bound states in the adiabatic 
approximation. We see that the symmetry present in Ref. |18j is partially broken 
since Eab = and only regions 2 and 3 are symmetric. The systems 87 Rb- 87 Rb- 6 Li 
(diamond) and 133 Cs- 133 Cs- 6 Li (circles) have three and four bound states, respectively. 
This is the experimentally relevant case [301 EH] > and our results demonstrates that 
also in a 2D setup, the 6 Li- 133 Cs system will have several bound states. 

It is known that the most favorable scenario for a spectrum with many bound 
states for two dimensional three-body systems is the symmetric energy case, where 
all the subsystems are equally bound, i.e. Eab — Eac = E B c [HI- However, 
this scenario seems less likely in current experiments. A very promising scenario 
is shown in figure [7J where the heavy-heavy system does not interact and the two 
heavy-light systems have the same binding energy. Except for the energy-symmetric 
case in Ref. [18], this is the most favorable scenario for a rich energy spectrum in 2D. 

Figure [8] shows scenario where the two-body energies are fully imbalanced. We 
have selected a generic case with energy scales related through Eab = IOEac and 
Ebc = O.IEac- This may appear difficult to set up experimentally, but fortunately 
this scenario does not seem to have any advantage over the others. In region 1, 
where our reference system has the heavy-heavy system more strongly bound than 
the others, the 87 Rb- 87 Rb- 6 Li and 133 Cs- 133 Cs- 6 Li systems have only one bound state 
each. In this energy configuration, the region 2 should be the most similar to region 
1 in figure [7J where the heavy-heavy system is not as bound as the others. However, 
regions 2 and 3 appear almost symmetric in figure [31 showing that both systems have 
two bound states each. Effectively, we see that the strongly bound heavy-light system 
changes the threshold of binding for the three-body system and effectively removes 
the most weakly bound three-body states. 
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Figure 8. Mass-diagram for an imbalanced energy system, defining the mass 
ratios thresholds for the existence of TV bound states for E^g = WEaq and 
EgQ = O.l-E^C- The numbers 1, 2 and 3 label three different regions according 
to figure [6] and the Roman numerals point out the number of bound states inside 
each solid curve. 



5. Summary and conclusion 

We consider three-body systems with zero-range interactions for general masses and 
interaction strength in two-dimensions (2D). The aim of this work is to study the 
limit where one particle is much lighter than the others (mc <C tua = tub) for non- 
interacting heavy particles, i.e. Eab = 0. This limit contains a rich energy spectrum 
and can be handled through the Born-Oppenheimer (BO) approximation. The BO 
approximation allows us to integrate out the light particle coordinate in order to derive 
the adiabatic potential between the heavy particles in the heavy-heavy-light system. 

The adiabatic potential is found to be the solution of a transcendental equation. 
We are able to find analytical forms for the asymptotic behavior for large and small R, 
where R is the distance between the two heavy particles. Comparing both asymptotic 
behavior and numerical solution of the adiabatic potential, we find that the analytic 
form precisely describes the adiabatic potential for any R and can be directly applied in 
2D three-body system calculations. The analytic adiabatic potential explicitly shows 
the mass-dependence. The analytic adiabatic potential becomes more attractive and 
suffers less screening as the mass of the light particle C decreases (mc — > 0). This 
potential explain the increasing number of bound states in this limit and we are able 
to demonstrate that the bound states accumulate both at R — > and R —> oo. We 
thus see an Efimov-like effect when mc = 0. However, we still caution that for any 
finite mc there are only a finite number of bound three-body states. 

The results obtained with the analytic adiabatic potential agree with numerical 
solutions of the momentum-space Schrodinger equation. However, approaching the 
three-body problem in 2D with the analytic form of the adiabatic potential has the 
huge advantage that it is much less computationally demanding and allows for easy 
access to rich regions of the mass-diagrams. Furthermore, we are able to show that 
the number of bound states calculated with the analytic adiabatic potential matches 
the estimate given by the semi-classical methods of the old quantum theory, i.e. the 
Bohr-Sommerfeld method. 

In recent work it was shown that the richest 2D energy spectrum is found when 
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the three-body mass-imbalanced system has all two-body subsystems bound with the 
same energy, namely Eab = E^c — Esc f° r m c "C tua ^ mg [T5]. However, as 
these conditions are difficult to obtain in experiment, we search for other configurations 
able to produce a rich energy spectrum. We numerically solve the momentum-space 
Schrodinger equation for general masses and interaction strengths in order to build 
mass-mass diagrams, which outline the expected number of bound state for a general 
three-body system. We find that systems with equal mass and non-interacting heavy 
particles, i.e. mc <C rriA = tub and Eab = 0',Eac = Ebc, are the most promising 
candidates to achieve a rich three-body bound state spectrum in 2D. Incidentally, this 
promising configuration is the one that can be handled with great precision in the BO 
approximation as we have shown. 

One of these promising configurations was recently reported as the experimental 
realization of 133 Cs- 6 Li systems [30l [31] . It even looks as if three-body bound states 
can be expected when the subsystem 133 Cs- 133 Cs is almost non-interacting. We can 
estimate the number of bound states for such a system in a two-dimensional setup 
through the mass-mass diagrams in this work. For 87 Rb- 87 Rb- 6 Li and 133 Cs- 133 Cs- 
6 Li we find respectively 3 and 4 bound states. It is very important to note that these 
numbers do not depend on the exact two-body energy in the 6 Li- 133 Cs subsystem. 
This two-body energy in the 2D setup are functions of the 3D low-energy scattering 
length [T] of the particular Feshbach resonance that is used in experiment to tune the 
interaction. However, as long as there is such a resonance, our results should apply 
when the system is squeezed into a two-dimensional geometry. 

For now, the system 133 Cs- 133 Cs- 6 Li prepared as in region 1 of figures [6] and [7] 
seems to be the most promising realistic combination to achieve a richer three-body 
energy spectrum in 2D. Some other mixtures under intense study at the moment with 
large mass imbalances are Lithium- Ytterbium [32] and Helium-Rubidium 33 . Of 
particular interest is that there are different isotopes available in those systems so that 
the quantum statistics of the components can also be changed. In the case of layered 
systems with long-range dipolar interactions as shown in figure Q] is also extremely 
promising for finding bound states [M] |3SJ 1311 12H I3H1 HOI 132], and we expect that 
some of these will have a universal low-energy character. The possibility of tunning 
the binding energy of each pair and performing experiments mixing molecules and 
atoms should open new avenues for even richer two-dimensional three-body spectra. 
Another interesting discussion is what the trace of the three-body parameter is in 
the corresponding 2D bound states. The three-body parameter determines the overall 
scale of the three-body spectrum in 3D and has recently been shown to be universally 
connected to the two-body van der Waals interaction [HI [42] [43[ SU [45l 06] [47] . Even 
though the 2D three-body problem does not require a three-body parameter, the fact 
that it seems to have a clear physical origin indicates that it could still play a role in 
2D. 

In order to experimentally observe the presence of these three-body bound states 
one should be able to use similar techniques to those used for the study of Efimov 
states in 3D, i.e. loss measurements [48j |3] and photo-association [49] [50]. It may 
also be possible to use RF techniques as in recent experiments studied two-body 
bound states and many-body pairing in two-dimensional Fermi gases [51] [29] . Another 
possible experimental signature of 2D three-body systems is through the momentum 
distribution and the two- and three-body contact parameters which appears as 
coefficients [52] [53] [54]. These coefficients depends sensitively on the presence of 
bound two- and three-body states. 
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